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The electron affinities e A(S) are measured for the two isotopes 32 S and 34 S (16752.9753(41) and 
16752.9776(85) cm -1 , respectively). The isotope shift in the electron affinity is found to be positive, 
e A( 3i S) - e A( 32 S) = +0.0023(70) cm" 1 , but the uncertainty allows for the possibility that it may 
be either "normal" [ 14( 34 S) > e A( 32 S) ] or "anomalous" [ 14( 34 S) < e A( 32 S) }. The isotope shift 
, is estimated theoretically using elaborate correlation models, monitoring the electron affinity and 

the mass polarization term expectation value. The theoretical analysis predicts a very large specific 
P_l ' mass shift that counterbalances the normal mass shift and produces an anomalous isotope shift, 

*vl( 34 S) — e A( 32 S) = —0.0053(24) cm -1 . The observed and theoretical residual isotope shifts agree 
with each other within the estimated uncertainties. 
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PACS numbers: 31.30.Gs, 32.80.Gc, 31.15.ve, 32.10Hq 



I. INTRODUCTION 

Photodetachment microscopy, which is the analysis of the electron interference pattern naturally produced when 
• i— i 1 photodetachment occurs in the presence of an electric field [J 0], was applied to a beam of 32 S~ ions and allowed 
to measure the detachment thresholds corresponding to different fine-structure levels of the negative ion S~ and 
(— i ■ the neutral atom S [3]. The electron affinity of Sulfur, at 2.077 eV, is well suited for detachment by a tunable 
O ' dye laser, which provides a third way of measuring neutral S fine structure, besides VUV spectroscopy of S I lines 
and direct fine-structure resonance spectroscopy. Dye laser photodetachment of S~ was also used as a probe of 
microwave- induced transitions of hyperfine-split Zeeman transitions, which lead to a measurement of the hypcrfinc 
structure of 33 S~ The fine structures of S~ and neutral S, with the definition of the six "fine-structure detachment 
■ thresholds", labelled A, B, C, D, E and F (in the order of increasing excitation energy) are displayed in Figure [TJ 
^\ \ In the present work, photodetachment microscopy is used to measure the electron affinity e A (threshold C) of the 
CNl ■ even isotopes 32 and 34 of Sulfur. The sensitivity of the method made it possible to record a significant number of 
34 S detachment events, even though Sulfur was produced from a chemical compound with no isotopic enrichment 
. and detachment occurred at very low energies above threshold (namely in the sub-meV range). The accuracy of 
photodetachment-microscopy based electron affinity measurements makes it possible to get an estimate of the isotope 
. shift, precise enough to make comparison with theory significant. 
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On the theoretical side, the ab initio calculation of electron affinities is a challenge. Various methods for evaluat- 
ing electron binding energies and affinities are discussed by Lindgren [5j], presenting five different techniques, from 
the Koopmans-theorem method, up to the density-functional-theory. Lindgren's survey focuses on the many-body 
perturbation approach and only mentions, without further discussion, a sixth method that can always be used: 
separate many-body calculations for the initial and final states, using some elaborate variational technique, like multi- 
configuration or configuration interaction methods. This approach can indeed be successful for small systems (see for 
instance ref. @ for 7 Li) and is precisely the one attempted in the present work, although the correlation balance is 
more difficult to achieve for a large number of electrons. 

The electron affinity of Sulfur has been estimated, after the pioneer work of Clementi et al. [7[, by Woon and 
Dunning who treated the second row atoms through multireference single and double excitation configuration 
interaction calculations and by Gutsev et al. [sj using the coupled-cluster method. In a benchmark study of ab initio 
and density-functional calculations of electron affinities covering the first- and second-row atoms, de Oliveira et al. fic| 
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concluded that the best ab initio results agree, on average, to better than 0.001 eV with the most recent experimental 
results. For heavy systems, the relativistic effects become crucial [TJ E3 and an accuracy better than 0.04 eV (10%) 
was difficult to achieve for the electron affinity of lead in the full relativistic approach [13| . 

The search for a possible variation of the fine structure constant a has renewed interest in developing reliable ab 
initio computational methods for atomic spectra (l4l . [l5j j . Theory versus experiment comparisons of atomic isotope 
shifts can serve as sensitive tests of our computational ability for some important electronic factors. Stimulated 
by photodetachment experiments [l6l . [l7| , theoretical calculations on isotope shifts in electron affinities have been 
attempted for Oxygen [l8|,[lj|, using the numerical MCHF approach. Beryllium was another interesting target [2(| but 
ab initio calculations remain scarce, requiring the calculation of two properties - the electron affinity and its isotope 
shift - both highly sensititive to correlation effects. For these systems, the limited population and the restricted active 
space concepts were used to build the configuration spaces. Adopting similar correlation models for the S/S - system 
becomes prohibitive and efficient reduction strategies need to be found. 

In the present work, various single- and double-multireference valence expansions sets are explored for a full 
non-relativistic variational optimization of the wave function through the MCHF procedure. Core excitations from 
these multircference sets are included afterward through configuration interaction. Computational strategies are 
investigated and developed by monitoring the electron affinity and the mass polarization expectation values difference 
between the neutral atom and the negative ion. 

The experimental work and theoretical calculations are presented in section [TTI and HTT1 respectively. The comparison 
of observation and theory is discussed in section HVl 

II. EXPERIMENTAL MEASUREMENT OF THE ISOTOPE SHIFT 

A. Experimental setup 

1. Ion beam and isotope selection 

Photodetachment microscopy was performed on a beam of S - produced by a hot cathode discharge in a mixture of 
98% Ar and 2% CS2. This commercial mixture had no isotopic enrichment, so the isotopes of Sulfur were produced 
with their natural abundance, which meant only 4% [2l[ of 34 S~, i.e. only a few pA. With such a barely measurable 
ion current, all the electrostatic settings of the ion beam and the alignment of the laser in the interaction region had 
to be done with our Wien velocity filter set on mass 32. Then the electric field applied in the Wien filter was reduced 
by the factor 4/VT7 appropriate for shifting from mass 32 to mass 34. We checked that the electron interferograms 
obtained with this setting were not a residue of 32 S~ signal on the wing of the maximum of mass-32 transmittance. A 
primary observation was that the peaks of the mass spectrum, when recorded on the total ion current signal (including 
a very visible contribution of SH~ at mass 33) appeared well separated, which meant a mass resolution of 70 at least. 
A more quantitative limit of the possible admixture of the 32 signal at mass 34 was given by setting the velocity filter 
at mass 33, and observing that this actually let no photoelectron signal emerge from the background of the electron 
image, even though we recorded the impacts on the detector for a time longer than the one needed to reconstruct 
visible interference rings at mass 34. One mass unit away from its maximum, the 32 S~ current was thus much lower 
than the 34 S~ maximum current. Since isotope 32 is 22 times more abundant than isotope 34, this means a factor of 
attenuation of 100 at least. Two units away from the 32 maximum, the attenuation must be even more complete, so 
the few 32 S _ detachment events remaining at mass 34 are certainly negligible with respect to the electron background 
due to parasitic collisions of the ion beam with the diaphragms or the residual gas and to the 34 signal itself. 

2. Laser photodetachment 

As in previous photodetachment microscopy experiments on Sulfur 0, l22j , laser excitation was provided by a CW 
ring dye laser (Spectra-Physics 380A) operating with Rhodamine 590 in 5% methanol and 95% ethylene-glycol. Single- 
mode operation was achieved by means of a pair of intracavity Fabry-Perot etalons. The cavity length is servo-locked 
by means of an external sigmameter (23[. Thanks to stabilization of the sigmametcr itself on the wavelength of a 
dual-polarization stabilized He-Ne laser, the frequency of the tunable laser can remain stable within a few MHz for the 
typically 20 minutes needed to record every photoelectron interferogram. The wave-number of the laser is measured 
by an Angstrom WS-U lambdameter, with an accuracy better than 10 -3 cm -1 . 
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B. Experimental data 

1. Photodetachment interferograms 

Figure [2] gives an example of a pair of interferograms obtained from a double pass of the laser beam on the ion 
beam. This double pass makes it possible to obtain Doppler-free measurements, by averaging the responses of both 
spots, for they correspond to symmetric residual deviations from 90° of the laser- ion intersection angle (c/. [24]). Due 
to the rarity of 34 S, even after an accumulation time of 2000 s, the number of electrons counted per pixel is 8 at 
its maximum. Nevertheless this is enough for the fitting program to find the centre and contour of each spot, and 
to calculate a histogram of the average number of electrons counted per pixel at a given distance from the centre. 
The obtained radial profile, as shown on Figure |3l gives a much less noisy picture of the interference pattern and 
shows its excellent correspondence with theory (even though the actual fitting procedure is done on the 2D electron 
distribution). The phase, i.e. 2tt times the number of oscillations in the interferogram, is the essential parameter for 
determining the initial kinetic energy of the electron. The interferometric accuracy so obtained, of the order of a few 
10" 3 cm -1 , is orders of magnitude better than what a measurement of the spot diameter would provide. The latter 
would actually pertain to the domain of classical electron spectrometry, the accuracy of which is seldom better than 
1 meV, or a few cm . 

2. Data analysis 

In principle, photodetachment microscopy does not require series of photodetachment images to be recorded to get 
a measure of the electron affinity e A. Calculating the difference between the photon energy and the measured electron 
kinetic energy in a single experiment would give the result. However, any discrepancy between the expected and 
actual values of the electric field in the photodetachment region is able to produce a systematic shift of the measured 
photoelectron energies. This being a constant relative error, extrapolating the measured electron affinity down to zero 
initial kinetic energy, i.e. to the detachment threshold, provides a method for avoiding the electric field uncertainty 
[13, HBj]. Series of measured electron affinities obtained both for isotope 32 and isotope 34 are represented on Figure 111 

Extrapolation of the measured e A values down to zero being the leading idea, we had to admit that the few 
experimental points obtained would not be enough to determine the slope of the linear regression, for the 34 S case, 
with a satisfying accuracy. The idea was thus to make 32 S and 34 S measurements in similar experimental conditions, 
and make the linear regression with a constraint of similarity on the slopes, in order to set the 34 S slope with improved 
accuracy. The experimental data shown on Figure2]for 32 S were actually taken during the same runs as the S ones. 
The difference between the 34 S and 32 S slopes may thus be constrained by a normal distribution with a characteristic 
width of 0.5%, which is an estimate (on the larger side) of the typical slope variations observed in past experiments 
done in similar conditions. As a matter of fact, fitting the data with this constraint yields nearly identical slopes of 
0.18 % and 0.16 % for 32 S and 34 S respectively, the visual consequence of this result being that the two regression 
lines drawn on Figure [4] appear nearly parallel. 

3. Experimental results 

The electron affinities e A(S), which are the ordinates at zero energy of the two lines drawn on Figure [4] are 
16752.9753(41) and 16752.9776(85) cm" 1 for 32 S and 34 S respectively. The error bars given here take the statis- 
tical distribution of the data into account at a 2a level together with a possible ±10~ 3 cm" 1 systematic error on 
wavenumber measurements. The value of e A( 32 S) incorporates all our previous work that led to the published result 
of 16752.9760(42) cm" 1 [3]. The subsequent discovery that transverse magnetic field effects were actually negligible 
[2r| was applied to these former data, but accounts only for —0.0001 of the —0.0007 revision of the most probable 
value of e A( 32 S) down to 16752.9753(41) cm" 1 ; e A( 32 S) remains the most accurately known of all electron affinities. 

The isotope shift e A( 34 S)- e A( 32 S) is found to be +0.0023(70) cm" 1 . Having a better accuracy on the difference than 
on the least-well known of both electron affinities is a logical consequence of the experimental method. The direct 
comparison of the apparent electron affinities at similar detachment energies naturally provides a good accuracy on 
the isotope shift, even though the linear regression to the actual value of e A suffers from additional experimental 
unknowns. The strong covariance of the electron affinities is reinforced by the inclusion of the possible systematic 
error on wavelength measurements, which is, by definition, the same in both cases. Numerically, the final correlation 
found between the obtained variances of the electron affinities is +0.57. 
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III. THEORY AND AB INITIO CALCULATIONS 



A. Theoretical Isotope Shift 



Adopting the (A' > A) convention where A is the mass number, the isotope shift on the electron affinity defined as 

IS(A', A) = S e A= e A(A') - e A(A) (1) 

is expressed as the sum of the normal mass shift (NMS), specific mass shift (SMS) and field shift (FS) contributions 

IS(A', A) = S e A NMS + 6 e A sus + S e A FS . (2) 

Introducing M for the nuclear mass and X for the chemical element, the two first terms that constitute the mass shift 
can be written in atomic units (m e = 1 and h = 1) as 
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This expression is correct to first order in n/M, where /i = m e M / (m e + M) is the reduced mass of the electron with 
respect to the nucleus. For a positive AS SMS difference, the NMS and SMS interfer negatively due to the relative 
signs of the mass factors in equation ([3]) [20] . It is easy to show that the degree of cancellation between NMS and 
SMS is basically governed by the mass-independent difference [ e A(oo) — ^-AS am ] . The two atomic masses of 32 S 

(31.972 071 00 u) and 34 S (33.967 866 90 u), taken from the AME2003 compilation of Audi et d. [H are converted 
into nuclear masses by subtracting the electron mass contribution (0.027 %)[40|. 
The field shift (FS) can be estimated from 



S e A FS = (hc)4iT p(0)* R - p{0) 
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where p(0) jv_r is the non-relativistic spin-less total electron density 
/(Z) 32 " 34 = 0.014823 cm" 1 / 1 " 12 is taken from Aufmuth et d. 
relativistic density for a point nucleus. The (r 2 )^ 2 = 3.2608(18) fm and 
from Angeli [3(3]. 



(r) (in a 3 ) calculated at r = [HJ]. The factor 
and corrects for the fact that we use the non- 



r /34 — 



3.2845(21) fm values are taken 



Note from eq. (|6]) that with a positive variation of the rms nuclear charge radii, i.e. S( 



„2\AA' — 



'r 2 )A' 



> 0, 



the FS has the same sign as the NMS if and only if the electron detachment [X~ —> X) is accompagnied by an increase 
of the electron density at the nucleus (Ap(0) > 0). For a system like Sulfur, the FS is expected to be much smaller 
than the mass shift. Therefore, our computational strategy is dictated by the description of the electron affinity and 
ASsms, although the FS is taken into account in the present analysis (see section HV)) . 



B. Computational method 

We use the numerical multi-configuration Hartree-Fock approach (MCHF) describing the atomic wave function as 

* = Y,CiHliLS) (7) 

i 

where {^(-jiLS)} is an orthonormal set of configuration state functions (CSF) that are symmetry adapted linear 
combinations of Slater determinants [3l|,[32j]. In this method, the radial functions {P n i(r)} defining the orbital active 
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set and the mixing coefficients {ci} are variational. The configuration interaction (CI) method solves the eigenvalue 
problem in a CSF basis built with a fixed preoptimized orbital set. 

For any differential property that is estimated from the difference between two calculated diagonal properties using 
the variational approach, much care must be taken to obtain a good balance between the two states. It becomes even 
more difficult when the latter belong to systems with different numbers of electrons. In such situations, the ATSp2k 
package [33| is an efficient tool thanks to its flexibility. In particular, the fully implemented limited population (LP) 
and multi-reference (MR) [34l |35| approaches for building the configuration space offer systematic ways of including 
and monitoring correlation. While the LP configuration space is built by allowing single (S)-, double (D)-, triple (T)- 
and possibly higher excitations, from a single reference configuration state function, restricted by orbital occupation, 
the MR method is usually limited to SD excitations from a larger set of configuration states. The LP method has 
been successfully used to calculate the isotope shifts in the Oxygen electron affinity [H, [l9[ . Although both the LP 
and MR correlation models have been explored, the present theoretical discussion is limited to the multireference 
calculations. 



C. The experimental electron affinity as a guideline 

In our approach, the experimental electron affinity is used as a guideline to set efficient pathways in the variational 
configuration spaces. Our non-relativistic approach targets electron correlation. In this context, it is useful to get a 
reference non-relativistic e A value. Introducing the observed average energy levels 

T _ Ej(2J+ l)Ej 

for both S 3p 4 3 P and S" 3p 5 2 P°, the average experimental electron affinity that would be measured if not resolving 
the fine structure thresholds, is estimated from 

*AZ = E( S 3 P) _ E(S- 2 n = ^ +3 ^ + 5 ^ - (2 ^ 2 ; 4 ^ 2) . (9) 

This average energy can be rewritten in terms of some deviation to the observed electron affinity e A exp = (E 2 — £"3/2) 
(arrow C in Figure [TJ 



e A AV _ e A , SjE, - E 2 ) + (E ~ E 2 ) (E 1/2 - E 3/2 ) 



Using e A exp = 16752.9753(41) cm 1 reported in section Hi B 31 for 32 S and the observed fine structures of S and S 

the Sulfur average experimental electron affinity e Af x v p = 16787.55 cm" 1 = 2.081 391 eV = 0.076 489 702 E h is 

obtained from eq. (fTU|) . 

When adopting a Breit-Pauli description [32], [36[ of atomic structures, the total binding energy E(LSJ) of a level 
is expressed in first order perturbation theory as 

E(LSJ) = E™ + E» F +E F SJ , (11) 

i.e. as the summation of the non-relativistic total energy {E ), the non-fine structure relativistic shift (E ) and 
the J-dependent fine-structure correction (E F ). Since the fine structures of both S 3p 4 3 P and S~ 3p 5 2 P° are 
washed out by the averaging process (fTU|) to get e A^, a reference non-relativistic electron-affinity e A^ is estimated 
by subtracting the corresponding theoretical non-fine structure contribution AE NF = E NF (S) — E NF (S~), from the 
above experimental average electron affinity 

*A%}= e Af^ — AE NF . (12) 

The non-fine structure contribution calculated in the single configuration Hartree-Fock approximation, AE^ F = 
—5.36251 10~ 4 E^, produces a non-relativistic electron- affinity value of e A^Jj = 0.077 026 E^. Note that the latter 
value is in line with the estimation of the "non-relativistic experimental" electron- affinity (0.076 939 Eu) calculated 
from the scalar contribution of de Oliveira et al. [Toj , who found an excellent general agreement for electron affinities 
of first- and second-row atoms. 
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D. The Multi-Reference approach 

1. MR-MCHF calculations 
In the multi-reference approach, one first defines a zeroth-order set of CSFs labelled MR 

MR = {$i( 7 iLS7t), $ 2 ( 72j L5^), . . . , <S> m ( lm LSn)} (13) 

that includes the dominant interacting terms in the description of a given atomic state. This zeroth-order set is 
then expanded to capture major correlation effects. Useful expansions are built by allowing all single and double 
excitations from a multi-reference (MR-SD) set within a given orbital active space. From a practical point of view, 
these expansions are generated using LSGEN [37J that produces the desired list of configurations, containing the 
complete set of CSFs for a given LSir - total symmetry. 

A good valence correlation MR-SD expansion for S~ is based on the multireference set 

MR(S") = {1, 2} 10 {3s, 3p} 5 {3, 4} 2 . (14) 

The notation is inspired from the LP approach: the multireference set (Ti"4|) is composed of all CSFs having the required 
symmetry (here 2 P°), with ten electrons {I, 2} 10 forced to occupy the n = 1 and n = 2 shells (i.e. a ls 2 2s 2 2p 6 closed 
core). In the MR space, the seven valence electrons should describe the dominant configuration 3s 2 3p 5 , but they are 
also free to reorganize themselves in the n = 3 and n — 4 subshells with only one occupation constraint: a minimum 
of five electrons should be either 3s, or 3p, as explicitly stated through the notation {3s, 3p} 5 . But even with a 
closed core {1,2} , the computational effort is gigantic. Introducing the \n max lmax] notation for the orbital active 
set, the size of the expansion generated with six correlation layers (MR-SD \9k~\) reaches 1 895 416 CSFs. Moreover, 
such a strategy is not efficient, a large number of components being negligible. An interesting approach uses the 
"Multi-Reference Interacting" (MR-I) CSF-space defined as the union of the original set of CSFs that belong to the 
MR, and all CSFs that directly interact with at least one component of the MR, i. e. 

$,( 7 , LS) G MR-I & 3$ fe G MR with ($,( 7l LS) \H\ $ fc ( 7fc LS)) ^ , (15) 

for any one-electron radial function basis set {P n i(r)}. This selection constraint depends on the coupling ordering of 
the subshells. The conventional coupling hierarchy is a sequential one corresponding to the coupling of each subshell 
angular momenta to the previous intermediate coupling angular momenta, from left to right [38j . for the natural 
subshell ordering (n and I increasing). But this is not always the most efficient representation. It is indeed common 
knowledge that the most strongly interacting momenta should be coupled first to get the best physical picture of the 
resulting levels pattern. On this basis, we use the reverse order of orbitals, coupling sequentially the subshells by 
decreasing n and I. 

Like the MR-SD space, the MR-I configuration set includes at most double excitations with respect to the reference 
but the building rule (fl~5|) reduces drastically the size of the expansions in comparison with the MR-SD sets. For 
example, the list of 1 895 416 CSFs discussed above is reduced to 525 111 CSFs in the MR-I|~9fc] model. The 
corresponding configuration reduction is much smaller (18 576 — > 13 973) since this reduction only arises from the 
intermediate coupling constraints associated with (|15[) and not from additional orbital occupation number selection 
rules. From a practical point of view, the MR-SD expansions are reduced according to the building rule (TTBl to 
produce the desired MR-I lists using the lsreduce code integrated in the ATSP2K package [33l |. 

Similar to (TJH , we explore the MR set 

MRl(S) = {l,2} 10 {3s,3rf 4 {3} 2 , (16) 

for S. To discuss the delicate balance between the negative ion and the neutral atom, we introduce for the latter a 
second model, more correlated, based on the following MR set 

MR2(S) = {1, 2} 10 {3s, 3p} 4 {3} x {3, 4} 1 . (17) 

The orbital active sets are extended up to n = 9, limited to I < 7 (i.e. fc-orbitals). We use multireference sets with 
at most two excitations in the 3d subshell. Allowing more than one 3d electron is indeed necessary for all important 
intermediate couplings to appear in the set of CSFs satisfying dl~5l) . 

The MR set (|T4l contains 157 CSFs. Some sublists of this complete configuration space are selected according to 
their impact on the energy, mass polarization and density at the nucleus and used to check the consistency of our 
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TABLE I: Mean radius of spectroscopic orbitals in atomic units of length (ao) obtained for S and S with MR1 and MR2 
including 6 correlation layers (9k). 







S 3p 4 A P 




S" 3p b 


^ po 


nl 


HF 


MR1-I[9fcl 


MR2-If9fel 


HF 


MR-I|9fc] 


Is 


0.09715 


0.09715 


0.09715 


0.09715 


0.09715 


2s 


0.47577 


0.47245 


0.47248 


0.47585 


0.47276 


2p 


0.44104 


0.44061 


0.44061 


0.44106 


0.44037 


3s 


1.72072 


1.71439 


1.61071 


1.77672 


1.77863 


3p 


2.06072 


2.02942 


1.84593 


2.32369 


2.25355 



results. By investigating the impact of a seventh correlation layer (10k) using multireference subsets for S , this extra 
layer is estimated to contribute less than 5.10 E^ to both e A and ^-AS SM s. 

With the valence correlation expansions based on MRs (|14[) . fp~6|) and (fl~7|) . we choose to vary all orbitals in 
the MCHF approach, the frozen-core approximation being considered a priori artificial. The mean radii of the 
spectroscopic orbitals of S and S~ are reported in Table U and compared with the HF ones. As expected, we observe 
an overall stability of the n — 2 orbitals and a larger variation for the (n — 3) valence shells, while the Is orbital 
remains very similar in all calculations. Even if the (ls 2 2s 2 2p 6 ) core is kept closed in the MCHF expansions, some 
correlation orbitals extend into the inner region of the atom to improve the description of the total wave function. 
In the Sulfur MR2-I[9fc] calculation for example, the MCHF optimization involves 39 numerical correlation orbitals 
from which the resulting 6p, Id, 8s and 9/ functions can be qualified as "inner" orbitals by looking at their mean 
radius ((r) n i < ao). Although they still describe the (inner region) valence correlation, they lie in the correct region 
for estimating core- valence correlation effects through configuration interaction, as presented in the next subsection. 



2. Open-core configuration interaction calculations 

We add to the valence configuration lists, core-valence mono- and multi-reference SD expansions (MR-CV-SD) 
created by allowing at most one hole in the core but keeping the Is shell closed and inactive. Core- valence excitations 
generate much larger lists of configurations than equivalent valence expansions. 

For keeping the size of the expansions tractable we use the following procedure. First we sort the configurations of 
the original MR-I valence eigenvectors into decreasing order by their configuration weights. The latter is defined as 
the weighted contribution of the CSFs belonging to it: 




and are reported in Table [TJl for S~ and for the two valence models (MR1 and MR2) used for S. Following this 
hierarchy, we define p reference subsets MR P containing the first p configurations in the sorted list. Secondly, we build 
the corresponding MR p -CV-SD spaces and keep only the CSFs interacting directly with the complete MR. We denote 
unambiguously the open-core CI calculations - all performed with six correlation layers - MR-I/CVp. 

Table Hill reports the electron affinity theoretical values calculated with the MR (S _ ) and MR2 (S) models. The 
largest configuration interaction calculation remains feasible for Sulfur (MR2-I/CV15), but the computational limits 
are definitely exceeded in the negative ion (MR-I/CV31). For S~ indeed, a larger calculation than MR-I/CV14 
requires truncating the expansions [4l|. To construct the MR-I/CVp spaces for p = 20 and p — 31, we first omit the 
CSFs with a |ci| < 1.10 -6 in the preceding expansions (p = 14 and 20), with an impact smaller than 10 -6 E^ on both 

energy and ^-AS SM s- These lists are then completed by adding the CV expansions of configurations 15 — 20, and 
21 — 31, respectively. 

Table IIIII displays smooth convergence trends along its diagonal but the largest expansions used for both S and 
S~ definitely underestimate the non-relativistic experimental electron affinity ( e A^!3 = 0.077 026 E u ), indicating 
that the neutral system is too correlated with respect to the negative ion. We do not report the corresponding 
table for MR(S _ )/MR1(S) models (having 4 columns instead of 15) that displays a good convergence toward e A^Jj. 
For a given correlation model, the AS SMS parameter is calculated with the wave function expansions that bring the 
theoretical electron affinity value as close as possible to the e A^- reference value. This approach is supported by 
the strong correlation observed between the total energy and the S SM s parameter, as discussed below. When adding 
configurations, one by one, in the MR p of neutral S, we look for the corresponding model in S~ that gives the best 
energy balance. These values are underlined in Table HTT1 and only appear in the first three columns corresponding to 
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TABLE II: Weights (w) of the configurations composing MR1, MR2 of S and the MR of S~ in the corresponding MR-I[9fc] 
wave-functions. # denotes the configuration index. 



S 3p 4 3 P 








S" 3p l 


'I po 






MRl-ir9fcl 


MR2-I[9fc 


1 








MR-I 


'9k 








config. 


w 


config. 




w 


# 


config. 




w 


# 


config. 




w 


— ; — 77— ri 

3s 2 3p 4 


0.9567 


— ■ — rr— t 

3s 2 3p 4 





9151 


1 


3s 2 3p b 





9382 


16 


3p b 3d 2 





0296 


3s3p 4 3d 


0.1605 


3s 2 3p 3 4p 





2652 


2 


3s 2 3p 3 3d 2 





1773 


17 


3s 2 3p 3 3d4s 





0292 


3s 2 3p 2 3d 2 


0.1568 


3s3p 4 3d 





1573 


3 


3s 2 3p 3 4p 2 





1249 


18 


3s 2 3p 3 4s 2 





0291 


3p 4 3d 2 


0.0498 


3s 2 3p 2 3d 2 





1508 


4 


3s3p 5 3d 





1045 


19 


3p 5 3d4d 





0288 






3s3p 4 4s 





1256 


5 


3s 2 3p 4 4p 





1030 


20 


3s 2 3p 3 4s4d 





0287 






3s3p 3 3d4/ 





0655 


6 


3s3p 4 4s4p 





0748 


21 


3p 6 4p 





0252 






3s3p 3 3d4p 





0632 


7 


3s 2 3p 3 3d4d 





0664 


22 


3p 5 4p 2 





0239 






3s 2 3p 3 4/ 





0531 


8 


3s3p 4 3d4/ 





0647 


23 


3p 5 4d 2 





0203 






3p 4 3d 2 





0492 


9 


3s3p 4 3d4p 





0626 


24 


3p 5 4s 2 





0200 






3p 5 4p 





0310 


10 


3s3p 4 4p4d 





0599 


25 


3s 2 3p 3 4p4/ 





0173 






3s 2 3p 2 3d4d 





0247 


11 


3s 2 3p 4 4/ 





0537 


26 


3p 5 4/ 2 





0124 






3p 4 3d4s 





0142 


12 


3s3p 5 4d 





0447 


27 


3s3p 5 4s 





0074 






3s3p 4 4d 





0127 


13 


3s 2 3p 3 4d 2 





0439 


28 


3p 5 3d4s 





0021 






3p 4 3d4d 





0085 


14 


3s 2 3p 3 4/ 2 





0401 


29 


3p 5 4s4d 





0020 






3s 2 3p 2 3d4s 





0030 


15 


3s3p 4 4d4/ 





0322 


30 


3p 5 4p4/ 





0018 




















31 


3s3p 4 4s4/ 





0010 



TABLE III: Electron affinity ( e A, in E^) versus the number of configurations (p,p') in MR2 P 
absolute energy and total number of CSFs (NCSF) of each model MR-I/CVp is given in the first 
The configurations are taken in the order of increasing weight (see Table ITT)) . Underlined are 
agreement with e A' 



for S and MR p / for S". The 
lines and columns of the table, 
the values of e A in reasonable 



NR 



0.077 026 Ei. 



p \p 








S 


1 


2 


3 


4 


15 






NCSF 




235 971 


355 354 


537 163 


681 582 • ■ 


■ 2 407 805 


S" 


E(in E h ) 


-397.718278 


-397.722125 


-397.723707 


-397.724631 •■ 


• -397.726165 


1 




541 


780 


-397.794996 


0.076718 


0.072871 


0.071289 


0.070366 • ■ 


■ 0.068832 


3 




864 


954 


-397.797229 


0.078951 


0.075104 


0.073522 


0.072598 • ■ 


■ 0.071064 


4 




982 


233 


-397.797780 


0.079501 


0.075654 


0.074073 


0.073149 • • 


■ 0.071615 


5 


1 


088 


076 


-397.798944 


0.080666 


0.076819 


0.075237 


0.074313 • • 


■ 0.072779 


6 


1 


210 


344 


-397.799173 


0.080895 


0.077048 


0.075466 


0.074542 • ■ 


■ 0.073008 


10 


2 


623 


506 


-397.799929 


0.081651 


0.077804 


0.076222 


0.075298 • ■ 


■ 0.073764 


12 


2 


854 


430 


-397.800229 


0.081951 


0.078104 


0.076522 


0.075598 • ■ 


■ 0.074064 


14 


3 


175 


092 


-397.800372 


0.082094 


0.078247 


0.076665 


0.075741 • • 


■ 0.074207 


20° 


3 


839 


474 


-397.800532 


0.082254 


0.078407 


0.076825 


0.075901 • ■ 


• 0.074367 


3ja 4 


339 


910 


-397.800667 


0.082389 


0.078542 


0.076960 


0.076037 • ■ 


■ 0.074504 



"These lists are truncated (see text). The actual CSF numbers are 2 089 778 and 2 058 776 for p = 20 and 31, respectively. 



MR2-I/CVp (p = 1 — 3), all larger p > 4 values underestimating the electron affinity, even for the largest MR-I/CV31 
S~ calculation. The three associated S~ correlation models correspond to MR-I/CV1 (mono- reference), MR-I/CV6 
and MR-I/CV31, respectively. For the approach based on the less correlated model for S (MR1), we select, using the 
same criteria, the calculations MR-I/CV4, MR-I/CV6, MR-I/CV12 and MR-I/CV31 of S - for the four MRl-I/CVp 
(p = 1 — 4) of S, respectively. 



3. Valence and Core-valence results 



Table IIVI reports in two blocks the electron affinities and S SM s differences for the valence correlation models and 
for their open-core extensions, using the Sulfur MR1 and MR2 models, respectively. In the upper half of Table IIVI 
we compare the results of the valence correlation calculations using the MR1-I[9fc] model (see section UlID ip . with 
the values obtained from the four CI calculations based on the core-excited correlation models. In the MR(p,p') 
adopted notation, p refers to the model used for S while p' refers to S~. According to our approach, the theoretical 
electron affinity values are forced to align with the non-relativistic experimental value through the (p,p') selection, 
but there is no such constraint on AS SMS - Opening the core through the added CV expansion affects the AS SM s by 
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up to 12 percents. Since the MR1(4,31) calculation corresponds to the complete models, the extracted AS SM s is a 
priori reliable. We observed that the results are well aligned, for each system, when plotted in a total energy versus 
S SMS diagram. Furthermore, the relation is similar for S and S~. A close analysis of the convergence patterns of the 
MRl(p,p') results leads to a 10% uncertainty estimation on the calculated AS SM s- 

The second half of Table HVl displays the corresponding results using MR2 for Sulfur. A good consistency with the 
MR1 ASsms values is found for the valence and first open-core (MR2(1,1)) models, but the two larger core-excited CI 
calculations MR2(2,6) and MR2(3,31) bring unfortunate variations. The effect of the truncation at p = 3 of the Sulfur 
MR2-I/CV15 model is estimated by the MR2(4,31) configuration interaction calculation that slightly underestimates 
the NR electron affinity. One observes that this extension affects the AS SM s value by more than 15%. Furthermore, 
the complete calculation MR2(15,31) is unreliable, given its large underestimation of e A^. 

From all these observations, we reject the open-core CI models based on MR2. Indeed, the Sulfur model includes 
much more correlation than the one built for S~. If one goes from a complete model (MR2(15,31)) that strongly 
underestimates e A to a balanced model that adjusts the electron affinity, it must be through a too large truncation of 
the S expansions. The problem of underestimating the S _ correlation energy with respect to the neutral atom's one 
is then transferred onto another problem which is the lack of convergence for the latter system (S). 

The breakdown of the proposed open-core procedure using the MR2 model for S is probably due to the different 
nature of the total wave functions obtained for S~ and S. The MR1 approach produces more comparable wave 
functions for Sulfur and its negative ion, respecting the needed balance. Signs for a large difference between the 
MR1 and MR2 Sulfur wave function appear in the analysis of their representation, through the comparison of their 
respective spectroscopic orbital mean radii (see Table U) and of their configuration weights (see Table HI)) . 



TABLE IV: Number of CSFs (NCSF), total energy (E, in E h ) and S SMS parameters (in 
report the corresponding electron affinity ( e A, in E^) and AS SM s (in slq 2 ). The two sets 
multireferences MR1 and MR2 used for S (see text). For each set, the results from 
open-core CI calculations (MR(p,p')) are reported. 

S 3p 4 a P 



a^ ) for S and S~. The last two columns 
of results correspond to the zeroth-order 
the valence models (MR-I [~9fc] ) and the 



Model 



MRiTjWj 

MR1(1,4) 

MR1(2,6) 

MR1(3,12) 

MR1(4,31) 



MR2-I[9fc] 

MR2(1,1) 

MR2(2,6) 

MR2(3,31) 

MR2(4,31) 

MR2(15,31) 



NR exp." 



NCSF 



E 



43 276 
66 280 
181 851 
268 647 
408 152 



-397.673394 
■397.720746 
-397.722224 
-397.723184 
-397.723563 



-67.01749 
-66.64770 
-66.61536 
-66.60345 
-66.59150 



209 553 
235 971 
355 354 
537 163 
681 582 
2 407 805 



-397.674938 
-397.718278 
-397.722125 
-397.723707 
-397.724631 
-397.726165 



-67.01897 
-66.68918 
-66.65212 
■66.61143 
-66.59715 
-66.58303 



NCSF 



3p 
E 



T7 



po 



525 111 
982 233 

1 210 344 

2 854 430 
4 339 910 



397.751017 
■397.797780 
397.799173 
■397.800229 
■397.800667 



-67.10600 
■66.73163 
-66.71128 
-66.69588 
-66.69335 



525 111 
541 780 
1 210 344 
4 339 910 



397.751017 
-397.794996 
■397.799173 
397.800667 



-67.10600 
66.78570 
-66.71128 
-66.69335 



A(S 



-s-) 

AS SMS 



0.07762 
0.07703 
0.07695 



0.07705 



0.07710 



0.0885 
0.0839 
0.0959 
0.0924 
0.1018 



0.07608 
0.07672 



0.07705 



0.07696 



0.07604 
0.07450 



0.0870 
0.0965 
0.0592 
0.0819 
0.0962 
0.1103 



0.07703 



a non-relativistic electron affinity e A^R defined in section [III CI 



E. Theoretical fine structures 



The fine structure splittings are estimated by performing Breit-Pauli configuration interaction calculations, including 
the orbit-orbit interaction. The results are presented in Table [Vj At the Hartree-Fock level of approximation, a large 
discrepancy between theory and observation (~ 30 cm -1 ) is found for 2 P°/ 2 -3/2 °f ^ an d 3 Pi-2 of S. Exploring 
various models for building the zeroth-order non-relativistic wave function, we observe that the inclusion of term- 
mixing due to L5-breakdown does not improve the fine structure splittings. Valence correlation is definitely insufficient 
to get a satisfactory agreement, as reflected by the splittings reported in [39]. A "simple" correlation model - denoted 
SD in Table[V]-, based on single- and double-excitations up to 6g from a single configuration and allowing at most one 
hole in the 2p-subshell (i.e. keeping Is and 2s closed), improves significantly the agreement between the theoretical 
and observed fine structure splitting values. Unfortunately, this agreement is destroyed when progressively extending 
the reference space to our more elaborate correlation model, illustrating the difficulty of getting reliable ah initio fine 
structure splittings. 
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TABLE V: Comparison of Breit-Pauli fine structure splittings (cm 1 ) with observation. 



s- 


HF 


SD 


MR-I/CV31 


observed [3] 


2 no 
-"1/2-3/2 


-453.03 


-482.07 


-471.16 


-483.54 


s 


HF 


SD 


MR1-I/CV4 


observed [3J 


3 Pl-2 


-366.91 


-394.82 


-395.22 


-396.06 


3 Po-i 


-181.89 


-174.82 


-169.93 


-177.54 



IV. ISOTOPE SHIFT IN THE ELECTRON AFFINITY: COMPARISON OBSERVATION - THEORY 

The observed and theoretical isotope shifts in the Sulfur electron affinity, both determined in the present work, 
are compared in Table IVTl The observed isotope shift (IS) on the electron affinity of Sulfur is found to be positive, 
e A( 34 S)- e A( 32 S) = +0.0023(70) cm" 1 , but the uncertainty implies that it may be either "normal" [ e A( 34 S) > e A( 32 S)} 
or "anomalous" [ e A( 34 S) < e A( 32 S)]. The normal mass shift (NMS) is easy to estimate from the first term of eq. (|3|), 
NMS = 0.016898 cm -1 , using the observed electron affinity. The experimental value of the residual isotope shift (RIS) 
is obtained by substracting the NMS contribution to the total isotope shift, i.e. RIS = —0.0146(70) cm -1 . 

As far as theory is concerned, two sets of results are reported, omitting or including the core-valence excitations, 
as described in section UlI D 31 For each set, the electron affinity, the specific mass shift (SMS), the total mass shift 
(MS = NMS + SMS), the field shift (FS), the total isotope shift (IS = MS + FS) and the residual isotope shift 
(RIS = IS - NMS) are reported. The electron affinities are compared with the experimental non-relativistic electron 
affinity ( e A^) estimated as explained in section UlI CI and reported in the same table. The field shift is estimated 
from equation ([6]) by calculating the change in the electronic densities at the nucleus A/?(0)nr with the density 
program [281 ] - The error bars of the FS values arise from the uncertainty in the root mean squares of the nuclear 
charge distributions, converted in a 17% variation of S(r 2 ) AA . 

For the valence calculations, convergence with respect to the number of correlation layers is achieved in both MR1 
and MR2 correlation models. The results reported in the "valence" line are obtained by averaging the MR1-I|~9fc] 
and MR2-I [~9Af| electron affinities and AS S ms parameters reported in Tabic IIVI Their uncertainty is estimated as 
half the difference between the two averaged values. The residual isotope shift (RIS = —0.0191(3) cm" 1 ) compares 
satisfactorily with the experimental result (IS — NMS) = —0.0146(70) cm" 1 . 

Opening the core is a very difficult task in the MCHF procedure. The presented open-core CI results are limited 
to the MRl-based models, due to the breakdown of MR2. The 0.12% and 10% error bars reported on the electron 
affinity and the SMS values, respectively, are estimated from the convergence of the MRl(p,p') sequence of results 
(see Table HVl) . The resulting calculated residual isotope shift of —0.0222(24) cm" 1 , remains compatible with the 
experimental RIS and its error bars, but the theory-observation agreement is tenuous. Core excitations affect the 
SMS contribution by 15%, increasing the total isotope shift by a large factor (2.4). The field shift constitutes a small 
fraction (~ 2%) of the residual shift (SMS+FS) but constitutes an important contribution to the total isotope shift 

(is). 

Both theory and experiment agree with a strong cancellation between the specific mass shift (SMS) and the normal 
mass shift (NMS) contributions. Although the normal or anomalous character of the isotope shift in the Sulfur electron 
affinity cannot be strictly confirmed from the present work, the ab initio calculations are definitely in favour of an 
anomalous IS. One should keep in mind however that the theoretical error bars are estimated from an objective analysis 
of the correlation models but do not take into account core-correlation and relativistic effects that are systematically 
neglected. 
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TABLE VI: Experimental and theoretical electron affinity, total isotope shift (IS) and residual isotope shifts (RIS). For the ab 
initio calculations, the specific mass shift (SMS), the total mass shift (MS) and the field shift (FS) contributions are reported 
separately. All values in 10 _2 cm _1 . 

e A IS RIS 

observation 

exp. 1675298 +0.23(70) -1.46(70) 

NR exp. a 1690524 

theory 6 

SMS MS FS 

valence c 16867(169) 10+^ -1.94(2) -0.25(2) 0.036(6) -0.22(3) -1.91(3) 

+ core-valence d 16922(200) 10+ 2 -2.25(23) -0.56(23) 0.038(7) -0.53(24) -2.22(24) 

"non-relativistic electron affinity e Az?f defined in section [HI CI 
adopting the experimental NMS. 

Averaging the MR1-I["9fc] and MR2-I["9fc] results (see text). 
d using the MR1(4,31) results of Table HVl 
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FIG. 1: S~ fine structure detachment thresholds. Energies measured with wavenumbers in m — 1 . (figure taken from 
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FIG. 2: Double interferogram obtained from 34 S~ at a wavelength A=596. 89056(2) nm, in an electric field 291 V/m, for an 
accumulation time of 2000 s. The grey scale indicates the total number of electrons counted in each pixel. The data are 
recorded by means of a Quantar Technology Inc. particle detector of the series 3391 with the 2251 Image-Trak™ enhanced 
software. The presented image was reprocessed with Microcal™ Origin®. 




Radius (mm) 

FIG. 3: Average number of electrons counted per pixel, in the lower spot of Fig. [2] The continuous line is the radial profile 
calculated with the best-fitting parameters (as provided by an adjustment algorithm applied to the original 2D data). One 
among these parameters is the initial kinetic energy of the electron, here found to be 0.5819(37) cm' 1 (but the la error bar 
given here assumes no uncertainty at all on the electric field itself). 
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FIG. 4: Comparison of electron affinity measurements made in common series of experiments for 32 S (circles) and 34 S (squares), 
with the average trend for each isotope (dashed line is 32, continuous line is 34), assuming similar dependences of the apparent 
electron affinity as a function of the average electron kinetic energy. The data shown in Fig. [2] produce the experimental 34 S 
point at 0.509 cm -1 . 



